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Classical transition state theory (TST) is the cornerstone of reaction rate theory. It postulates 
a partition of phase space into reactant and product regions, which are separated by a dividing 
surface that reactive trajectories must cross. In order not to overestimate the reaction rate, the 
dynamics must be free of recrossings of the dividing surface. This no-recrossing rule is difficult (and 
sometimes impossible) to enforce, however, when a chemical reaction takes place in a fluctuating 
environment such as a liquid. High-accuracy approximations to the rate are well known when the 
solvent forces are treated using stochastic representations, though again, exact no-recrossing surfaces 
have not been available. To generalize the exact limit of TST to reactive systems driven by noise, 
we introduce a time-dependent dividing surface that is stochastically moving in phase space such 
that it is crossed once and only once by each transition path. 
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I. INTRODUCTION 

Transition state theory (TST) 1 : 2 : 3 : 4 : 5 : 6 : 7 plays a cen- 
tral role in chemistry because it provides both a tanta- 
lizingly simple approach to calculating chemical reaction 
rates and an intuitive interpretation of reaction mech- 
anisms through the identification of the bottleneck be- 
tween reactants and products. Although the rate could 
in principle be calculated exactly using known meth- 
ods, in practice, it is often out of reach because such 
computations are laborious, even if quantum effects are 
ignoredi&SiiS Therefore, the TST approximation to the 
rate has played a central role in kinetics, but its accuracy 
hinges on the optimal identification of the TS dividing 
surface. 11 Indeed TST is exact when the latter is crossed 
once and only once by each reactive trajectoryiSiiS While 
a constructive prescription for a no-recrossing TS di- 
viding surface has long been sought, it was only in 
the 70's that it was found in the special case of low- 
dimensional systemsiii*i2ii^ A general method that yields 
a no-recrossing TS dividing surface for reactions in arbi- 
trarily many (but finite) degrees of freedom has been de- 
scribed only recentlyAi 5 * 1 ^ 7 " Nonetheless, even approx- 
imate TS dividing surfaces located near the phase-space 
bottlenecks between reactants and products provide a 
portrait of the activated complex that has been used rou- 
tinely in the chemistry community. 

Reactions in solution are even more complex: the pres- 
ence of the solvent introduces a complicated many-body 
problem. Nevertheless, TST has been remarkably suc- 
cessful at providing accurate estimates of thermal reac- 
tion rates when a reasonable TS dividing surface can 
be identifiedjiSiiSiSSiSi Unfortunately, such estimates are 
seldom exact because the dividing surfaces available in 
the literature are not strictly free of recrossings. One 



approach lies in choosing the dividing surface so as to 
minimize the TST reaction ratei 12 : 22 : 23 i 24 i 25 The recross- 
ing problem is thereby alleviated, but usually not elim- 
inated. Alternatively, one can eliminate consideration 
of the nonreactive trajectories entirely by sampling the 
transition path ensemble (TPE) directly. 26 : 27 : 28 : 29 : 30 All 
such paths intersect the TS dividing surface and hence, 
as a collection, they provide a portrait of the activated 
complex. Using Monte Carlo approaches ^SLSI TPE can 
also be readily applied to complex molecular reactions in 
all-atom solvents but it does not provide the simple geo- 
metric picture of a dividing surface free of recrossings. 

The geometry of the finite dimensional no-recrossing 
TS dividing surface is described in this paper for chemical 
reactions that can be represented using stochastic mod- 
els of the solvent interactions. The moving TS dividing 
surface is most readily constructed in cases of uniform 
solvents such as those represented by the Langevin equa- 
tion as was summarized in a recent Letter i2I The present 
article expands that discussion, and extends the results to 
the colored noise case of the generalized Langevin equa- 
tion. It should be noted that a finite-dimensional-like TS 
dividing surface can be specified if an explicit infinite- 
dimensional model — viz. an infinite collection of har- 
monic oscillators — is used to represent the heat bath in 
the Langevin equationiSi 3 ^* 3 ^ In spite of its simplicity, 
this model leads to an excellent approximation to the 
rate constant 3 ^ 5 - that, remarkably, has been rederived 
without recourse to an explicit model of the heat bathi 3 ^ 
Within the finite-dimensional phase space of the solute, 
the reaction geometry has also been illustrated using the 
stochastic separatrix (i.e., the collection of all phase space 
points for which the reaction probability is 50%.) That 
work portends the present development^ 7 " 

In this paper, we introduce a generalized view of the 



2 



transition state that allows for the temporal variation of 
the solvent. While this has not yet been accomplished in 
an all-atom representation, we show here that if the sol- 
vent is modeled through the Langevin equation, a time- 
dependent TS dividing surface can be constructed that is 
strictly free of recrossings. This is achieved by identify- 
ing a privileged stochastic trajectory that remains in the 
vicinity of the barrier for all time. This Transition State 
Trajectory serves as a moving origin to which geometric 
structures in a noiseless phase space, such as invariant 
manifolds and a no-recrossing surfaceiiLii are attached. 
The time-dependent dividing surface thus obtained adds 
to the evolving understanding of the geometric structures 
separating reactants from products in the presence of 
noise. The picture presented here provides a detailed 
time-resolved description of the reaction dynamics and, 
in particular, gives a precise geometric interpretation of 
the collective reaction coordinate introduced in Ref. |3|J 
The outline of the present paper is as follows: Sec- 
tionlTTIrecapitulates the macroscopic model of the solvent 
dynamics represented by the Langevin equation. In sec- 
tion lllll we introduce the fundamental construction that 
allows us to separate the full dynamics into the noisy 
TS trajectory and the noiseless relative motion. The TS 
trajectory for a particle driven by a white noise bath is 
described explicitly in section IIVI The construction is 
simple and direct, and requires no more effort than that 
for a typical trajectory. Several statistical properties of 
the TS trajectory are also discussed in section llVl In ad- 
dition, the TS trajectory is used to to compute reaction 
probabilities and the stochastic separatrix. The general- 
ization of the TS trajectory to the chemically important 
case of colored noise is discussed in section Ivl 



II. PRELIMINARIES 

The influence of the solvent on a chemical motion or 
reaction coordinate, q, has been routinely represented us- 
ing the Langevin equation (LE) of motionj22iSi 



q a (t) = -V ? U{q a (t)) - Tq a (t) + &(t) , 



(1) 



or its generalization (GLE) to allow for colored 
noise, 38 i 39 i 40 i 41 i 42 i 43 i 44 i 45 In Eq.ffl the vector ^denotes a 
set of mass- weighted coordinates in a configuration space 
of arbitrary dimension N, U(q) the potential of mean 
force governing the reaction, T a symmetric positive- 
definite friction matrix and £ a {t) a stochastic force. The 
subscript a here and throughout the paper labels a par- 
ticular noise sequence £ a (t)- For any given a, there are 
infinitely many different trajectories q a (t) that are distin- 
guished by their initial conditions at some time to. Note 
that no restriction on the dimensionality of the configu- 
ration space is imposed and the derivations below work 
in one as well as in many dimensions. Although in Eq.^ 
the deterministic force is formally assumed to be derived 
from a potential U(q), it would be straightforward to 



include velocity-dependent forces from, e.g., a magnetic 
fields Meanwhile, the stochastic force £ Q (£) is assumed 
to be Gaussian distributed with zero mean. It is re- 
lated to the symmetric positive-definite friction matrix T 
through the fluctuation-dissipation theorem) 3 ^ 4 * 4 ^^ 4 ^ 



&(*)£(*')} =2k B TT6(t-t') 



(2) 



where the angular brackets denote the average over the 
instances a of the fluctuating force, and the Dirac 6 func- 
tion appears because white noise is local in time. 

The assumption of white noise in Eq. [5] is consistent 
with the LE: It requires that the dynamics of the heat 
bath, which determines the correlation time of the fluctu- 
ating force, takes place on much shorter time scales than 
the dynamics along the chosen system coordinates — viz., 
the reaction coordinate and any other solute or solvent 
modes coupled to it. This is not always the case in appli- 
cations relevant to chemistry where the dynamical time 
scales of the system and the bath are usually comparable. 
Nevertheless, the special case of white noise is treated in 
detail in the next two sections because it accommodates 
an explicit description of the relevant geometric struc- 
tures. The generalization of the moving TS structures to 
the GLE is formally straightforward though mathemati- 
cally more cumbersome. In the case of colored noise, the 
fluctuation dissipation theorem requires that the friction 
kernel be nonlocal in time, viz. contain memory^ 4 ^^ 
The TS trajectory is nevertheless computable and is pre- 
sented in section Perhaps surprisingly, each of the 
geometric structures associated with the GLE is analo- 
gous to that in the LE. 

A further assumption in this work requires that the 
the reactant and product regions in configuration space 
are separated by a potential barrier, as is the case for 
most chemical reactions. The position of the barrier is 
marked by a saddle point q\ = of the potential /7(g). 
In the absence of noise, the saddle point is a fixed point 
of the dynamics. The invariant manifolds that determine 
the relevant phase space geometrjiiSiii are attached to 
it. In the presence of a fluctuating force, the geometric 
structures fixed at the saddle point no longer determine 
the reaction. As shown below, however, there is a unique 
stochastic trajectory q\ (t) that can play the role of the 
saddle point and serve as the carrier of invariant mani- 
folds. 

Because the reaction rate is determined by the dynam- 
ics in a small neighborhood of the saddle point(^& the 
deterministic force in the LE of Eq. ^ can be linearized 
around the saddle point to yield 



q a (t) = Slq a (t) - Tq a (t) + ( a (t) , 



(3) 



where the first term in the RHS involves the symmetric 
force constant matrix, 



CI, 



d 2 U 



(4) 



q=q+ 
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Although any rectilinear coordinate system can be cho- 
sen, it is convenient to choose the one for which the force 
matrix, 



n = 



(u 



\ 



\ 



(5) 



is diagonal from the outset. In Eq. qi is the reaction 
coordinate, uj^ the barrier frequency and u>2,.. -Wn the 
frequencies of oscillations in the stable transverse normal 
modes q 2 , . . . ■ 



III. THE FUNDAMENTAL CONSTRUCTION 

The aim in the following is first to construct a no- 
recrossing surface for Eq. and then to generalize the 
construction to correlated noise. It is well known that 
in the absence of noise and of dissipation — which are 
closely related by Eq. EJ— the phase-space hypersurface 
(<7i = 0) is f ree °f recrossings4i£ii£»il. As shown be- 
low, the damping does not pose particular difficulties. 
The noise, however, does, because the fluctuating force 
leads the particle to move randomly back and forth, so 
that it will typically cross and recross any fixed divid- 
ing surface many times. We will overcome this difficulty 
by constructing a dividing surface that is itself randomly 
moving such as to avoid recrossing. To achieve this, the 
dividing surface must be constructed as a function of the 
fluctuating force. 

Given any two trajectories q a (t) and q\{t) under the 
influence of the same fluctuating force, we define the rel- 
ative coordinate, 



Aq(t) = q a (t) - qi(t) . 



(6) 



It describes the location of the trajectory q a (t) with re- 
spect to the moving origin q\ (t) . By Eq. [21 the relative 
coordinate Aq(t) satisfies the deterministic equation of 
motion 



Aq = nAq-TAq 



(7) 



and thus exhibits non-random noiseless dynamics. This 
result is indicated by the lack of a subscript a in Aq. Due 
to the presence of the damping in Eq. [7J the dynamics of 
the relative coordinate, though noiseless, is still dissipa- 
tive. Nevertheless, with respect to the relative dynamics, 
it is possible to specify a no-recrossing surface, as well as 
the invariant manifolds of the equilibrium point Aq = 0. 
These geometric objects in the noiseless phase space can 
be regarded as being attached to and propagating with 
the reference trajectory q\{t). They thus define moving 
invariant manifolds and a randomly moving no-recrossing 
surface in the phase space of the original, noisy system. 

The relative dynamics between any two trajectories, 
q a (t) and q\(t), is described by Eq. [7] Because it is 



always noiseless, a moving no-recrossing surface can be 
thought of as attached to any reference trajectory. How- 
ever, only a specific surface is relevant for the reaction 
dynamics: The crossing of the surface should indicate 
the transition of the trajectory q a (t) from the reactant 
to the product side of the barrier. It cannot serve that 
purpose if it is attached to a reference trajectory that is 
itself far away from the barrier. If the reference trajec- 
tory is chosen arbitrarily, it will typically descend into ei- 
ther the reactant or the product wells over time and thus 
lose its ability to carry a chemically meaningful dividing 
surface. Only a reference trajectory that remains in the 
vicinity of the barrier for all time without ever descend- 
ing on either side of it can carry a no-recrossing surface 
that describes the reaction in the same way that a static 
dividing surface in conventional TST does. Indeed, we 
will show below that for each instance of the noise there 
is a unique reference trajectory with this property. This 
is the Transition State (TS) Trajectory mentioned in the 
introduction. We will henceforth restrict the notation 
q\(t) to this particular privileged reference trajectory. 



The definition of the TS trajectory as "remaining close 
to the saddle point" for all time might at first appear 
somewhat vague. In the construction of the TS trajec- 
tory below, it will become clear that the exponential in- 
stability that is inherent in the noiseless dynamics as- 
sociated with Eq. [21 in both the distant past and the 
remote future must be absent from the TS trajectory. 
A precise meaning can be given to the notion in statis- 
tical terms as follows: Because it is stochastic, at any 
given time there is a probability distribution for the TS 
trajectory in phase space. This distribution is invariant 
under the time evolution given by Eq. [21 In mathemat- 
ical language, the TS trajectory represents an invariant 
measure of the noisy dynamical system^ That is, for 
any instance £ Q ( - ) of the noise, where the dot indicates a 
function of time, there is an instance q\ (t) of the 

TS trajectory that is given as a functional of the noise 
and is itself a function of time, t. It can be used to ob- 
tain the instance of the TS trajectory q\ (t + r) 

that will be found if the origin from which time is mea- 
sured is moved an increment t / 0. Alternatively, one 
can use the time-shifted noise £ Q (- + r) and compute the 
corresponding TS trajectory q\ £ Q (- + r) (t). The TS 

trajectory is uniquely characterized by the requirement 
that these two ways of computing the time shift agree: 

&[£,(•) (« + r)=g* [£,(•+ r)l (t). 



We have now given a general outline of the method. 
In the following sections, we will actually construct the 
TS trajectory, determine its statistical properties and de- 
scribe the invariant manifolds and the no-recrossing sur- 
face attached to it in detail. 
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IV. THE TRANSITION STATE IN WHITE 
NOISE 

Equation can be rewritten as a first-order equation 
of motion in the 2./V-dimensional phase space with the 
coordinates 



(8) 



where v = q, as 



:;,</: .ir,,m~(£° ) on 



with the 27V-dimensional constant matrix 

.i = i v i . (io) 



where I is the TV x N identity matrix. After diagonaliz- 
ing the matrix A, its eigenvalues are denoted by 6j and 
its corresponding eigenvectors by Vj . The Langevin equa- 
tion, Eq.|21 then decomposes into a set oi2N independent 
scalar equations of motion 



(11) 



where z a j are the components of z in the basis Vj of 
eigenvectors of A and ^ a j the corresponding components 
of (CU>)) T 



A. Relative dynamics 

The noiseless equation of motion for the relative coor- 
dinate Aq(t) in Eq. when lifted into phase space via 
Eq. |S]and expressed in diagonal coordinates, reads 



It is solved by 



Aij(t) =ejA Zj (t) . 



Azj(t)=e^Azj(0) 



(12) 



(13) 



with a set of integration constants Azj(0). Equation 1131 
leads to a solution of the form 



e-it 



C jV 3 e i 



for the relative coordinate vector Aq(t) , where vj are the 
configuration space parts of the phase space eigenvectors 
Vj and the Cj are arbitrary constants. To identify reactive 
trajectories, one must consider the behavior of Aq(t) in 
the remote future and in the distant past and then deter- 
mine those trajectories that pass across the barrier from 
the reactant to the product side. 

The simplest case occurs in one degree of freedom 
where ft reduces to the squared barrier frequency u>^ and 



r to a scalar damping coefficient 7. In this case, the ma- 
trix A can be diagonalized analytically. It possesses the 
eigenvalues 

^ = -\ (7- vW 4 ^) >° ( 15a ) 



and 



fs = 

2 



\ (7 + V / 7 2 + 4 ^b) 



< 



with the corresponding eigenvectors 



V,, = 



and 



V = 



(15b) 



(16) 



The phase portrait of the dynamics is shown in Fig. [Qi. 
The eigenvectors corresponding to the positive and neg- 
ative eigenvalues, respectively, span one-dimensional un- 
stable and stable manifolds of the saddle point. They act 
as separatrices between reactive and nonreactive trajec- 
tories. The knowledge of the invariant manifolds allows 
one to determine the ultimate fate of a specific trajectory 
from its initial conditions. It should be clear that line, 
Aq u = 0, in the quadrant of reactive trajectories acts as 
a surface of no recrossing. 

In the absence of damping (and in units where lu^ = 1), 
the invariant manifolds bisect the angles between the co- 
ordinate axes. The presence of damping destroys this 
symmetry. As the damping constant 7 increases, on ob- 
tains the limits, e s — > and e u — * 00. Therefore, the 
unstable manifold rotates toward the g u -axis, the sta- 
ble manifold toward the g u -axis. In the limit of infinite 
damping, the invariant manifolds coincide with the axis. 
Thus, the fraction of phase space on the reactant side 
that corresponds to reactive trajectories decreases with 
increasing friction. This is intuitively clear because if the 
dissipation is strong, a trajectory must start at a given 
initial position with a high velocity to overcome the fric- 
tion and cross the barrier. 

In multiple dimensions, the dynamics of the single 
unstable degree of freedom is the same as in the one- 
dimensional case shown in Fig. Transverse damped 
oscillations, as shown in Fig. QJa, must be added to this 
picture. Their presence manifests itself, for small damp- 
ing, through N — 1 complex conjugate pairs of eigenvalues 
(14) £j- For stronger damping, some of the transverse modes 
can become overdamped, as illustrated in figure so 
that further eigenvalues become negative real. In any 
case, there is exactly one unstable eigenvalue with pos- 
itive real part. This eigenvalue is actually real and cor- 
responds to the system sliding down the barrier. In all 
other directions, the dynamics is stable, and at least one 
eigenvalue is negative real. The dynamics in the distant 
past is determined by the eigenvalue or pair of eigenval- 
ues with the largest negative real part. 

The eigenvector corresponding to the most stable 
eigenvalue together with the unstable eigenvector span 
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FIG. 1: Phase portrait of the noiseless damped dynamics corresponding to the linear Langevin equation in (a) the unstable 
reactive degree of freedom, (b) a transverse stable degree of freedom, (c) an overdamped transverse degree of freedom. 



a plane in phase space in which the dynamics is given 
by the phase portrait of Fig. Because the dynamics 
in all other, "transverse" directions is stable, the sep- 
aratrices between reactive and nonrcactive trajectories 
that we identified for the one-dimensional dynamics to- 
gether with the transverse subspace form separatrices in 
the high-dimensional phase space. In a similar manner, a 
no-recrossing curve in the plane together with the trans- 
verse directions form a no-recrossing surface in the full 
phase space. 

If the most stable eigenvalues occur as a complex con- 
jugate pair, the leading dynamics in the remote past is 
given by an oscillation across the barrier that, because it 
is strongly damped, has a huge amplitude in the distant 
past and causes the system to cross the barrier infinitely 
often. This behavior is clearly unphysical. In fact, once 
the oscillations become large, nonlinearities of the reac- 
tion dynamics become relevant, and the ultimate fate of 
the trajectory can no longer be determined from the lin- 
ear approximation. This situation can only occur if the 
damping is strong and at the same time the undamped 
transverse frequencies are so large that the correspond- 
ing eigenmodes do not become overdamped. This rare 
situation is ignored below. 

An exception to the discussion above occurs if for any 
eigenvector Vj the component in the direction of the 
reaction coordinate is zero. In this case the excitation of 
the corresponding eigenmode has no impact on whether 
a trajectory is reactive or nonreactive. The unstable and 
the most stable eigenvalues must then be chosen among 
the eigenvalues with ^ 0, which we call the "rele- 
vant" eigenvalues. Because the unstable eigenmode cor- 
responds to the particle sliding down the barrier, it is 
always relevant. Similarly, there must be one relevant 
negative eigenvalue corresponding to the particle being 
shot at the barrier from infinity and coming to rest on 
it. The occurrence of irrelevant eigenvalues might seem a 
rather unlikely exception, but it does arise in the impor- 
tant special case of isotropic friction, where the friction 
coefficient r = jl is a scalar multiple of the identity 



matrix. In this situation, the friction does not couple 
different normal modes of the undamped dynamics, and 
all but two eigenvalues are irrelevant. In particular, our 
construction of a no-recrossing surface can always be car- 
ried out for isotropic friction. 



B. The Transition State Trajectory 

The influence of noise in the equation of motion, 
Eq. ^] can be accounted for by means of the retarded 
Green's function 



G let (t) 



if t < 
e^* if t > 



(17) 



which satisfies 

G Tet (t)-e j G xe \t)=S(t) . (18) 
It gives the general solution of Eq. ^2 as 



z aj (t) = c aj e^ + / G lct (t - T )£ aj (t) dr 

J —OQ 

= c aj e^+ f e e ^- T H aj {T)dr (19) 

J — CO 

with arbitrary constants c a j. The integral 



4,-(*) = 



ej (*-"r) 



I e^ T £ aj (t + T)dT 

J — OC 



(20) 



exists with probability one if Re ej < 0, i.e., for the stable 
components. Only for those, therefore, is the solution, 
Eq. ^jj meaningful. 

The first term in Eq. 1191 tends to zero as t — > oo, but 
it diverges exponentially as t — > —oo. The only solu- 
tion to Eq. II II that remains bounded in the remote past 
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is obtained by setting c a j — 0. But the TS Trajectory 
has been defined as the unique trajectory that remains in 
the neighborhood of the saddle point at all times. It now 
becomes clear that this condition indeed determines the 
stable components of the TS trajectory uniquely and that 
they are given by Eq. [201 Because they depend on the 
noise £ a j(-), the z„j(t) are stochastic functions of time. 
From the second form given in Eq.EOI it should be obvi- 
ous that the probability distribution of z^j is stationary 

in time if the distribution of £<*(•) is itself stationary. 

The solution in Eq.^5]fails f° r the unstable component 
of the trajectory, where Ree^ > 0. To compute that 
component, we use the advanced Green's function 



G adv (t) 



to obtain the solution 



if t < 
if t > 



zj(t) =c aj e^ +z t aj (t) 



(21) 



(22) 



with constants c a j and 



= e-^Z aj {t + T)dT 

Jo 



(23) 



The position in Eq. [221 is bounded in the remote future 
only if c a j = 0. Therefore, z a ,(t) given in Eq. [221 can be 
identified as the component of the TS Trajectory in the 
unstable direction. Notice that the stable components in 
Eq. [201 at time t depend on the past of the fluctuating 
force, whereas the unstable component in Ea. l23l depends 
on its future. 

Through Eqs. [201 and [23 all components of the TS 
trajectory zj(i) in phase space have been specified. 
The TS trajectory q\{t) in configuration space can 
be obtained because the transformation connecting the 
position-velocity coordinates and diagonal coordinates is 
known. The consistency of the calculation can then be 
checked by verifying that the velocity obtained from the 
coordinate transformation is actually the time derivative 
of the position coordinate. 

A specific instance of the TS trajectory for a system 
with two degrees of freedom (and a four-dimensional 
phase space) is shown in Fig. [21 It can be observed how 
the trajectory randomly wanders around in the vicinity 
of the saddle point without ever leaving it. The veloc- 
ity space projection of the TS trajectory appears much 
more noisy than its configuration space projection. This 
results because the environment is assumed to cause a 
fluctuating force in Eq. Q that couples directly only to 
the velocity, but not to the position. It should be clear 
from the above discussion, however, that this condition 
can easily be relaxed. In the absence of noise and of 
damping, any trajectory would exhibit an elliptical orbit 
in the projection onto the stable degree of freedom This 



behavior is still visible in the noisy situation of figure 
where a slow change of energy, corresponding to the size 
of the ellipse, is superimposed over the elliptic motion. 

As described in section III, the TS trajectory serves 
as the origin of a moving coordinate system, and the ge- 
ometric structures in the noiseless phase space that we 
have described in section llV Al are attached to it. This 
is illustrated in Fig. [01 for two degrees of freedom. This 
case already exhibits all salient features of the dynam- 
ics in arbitrary high dimensions. The figure compares a 
specific instance of the TS trajectory with a transition 
path under the influence of the same noise. It can be 
seen how the transition path approaches the TS trajec- 
tory from the reactant side, remains in its vicinity for a 
while and then wanders off to the product side. Because 
the moving invariant manifolds of the TS trajectory are 
known, it can be predicted already at time t = that 
the trajectory actually will lead to a reaction instead of 
returning to the reagent side of the saddle. From both 
the time-dependent plots and the time-independent pro- 
jections in Fig. [OJi and [0t it is clear that the transition 
path crosses the space-fixed TS surface q u — = sev- 
eral times before it finally descends on the product side. 
The moving TS surface Aq u = 0, by contrast, is crossed 
only once, at the reaction time < rca ct = 8.500. That this 
is actually the case can be seen from the curves in the top 
faces of each column, which indicate the noiseless rela- 
tive motion Aq(t), Av(t) between the TS trajectory and 
the transition path. In the relative coordinates we find, 
as expected, the hyperbolic motion known from Fig.n m 
the unstable degree of freedom, a damped oscillation in 
the stable transverse degree of freedom, and a superpo- 
sition of the two in the configuration space projection in 
Fig. El:. 



C. Statistical properties of the TS trajectory 

If the noise £ a (t) is assumed to be Gaussian with zero 
mean, it is clear from Eqs. [201 arj d [201 that the same is 
true for all components z a a of the TS trajectory and, by 
extension, also for its position and velocity coordinates. 
To specify their joint probability distribution completely, 
it thus remains to calculate their time autocorrelation 
and cross-correlation functions. 

The most general fluctuating force £ Q (t) satisfying the 
fluctuation-dissipation theorem, Eq.[21 can be written as 



&»(*)= ffW tt (f) 



(24) 



in terms of M realizations of standard white noise 
W a (t) = (W a i(t), W aM (t)) T normalized to 

(W ai (t)W aj (s)) a = S(t - s) % (25) 

and anJVx M matrix g of coupling strengths satisfying 

gg T = 2k B Tr. (26) 
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FIG. 2: A random instance of the TS trajectory in a system with N — 2 degrees of freedom, projected onto (a) configuration 
space, (b) velocity space, (c) the reactive degree of freedom, and (d) the transverse degree of freedom. Units have been chosen 
so that ojb = 1, IcbT = 1. The transverse frequency is ui2 = 2, and friction is isotropic, T = yl with 7 = 0.5. 



The statistical properties of the noise are completely 
specified by Eq. If M > N, a representation can 
be found that has fewer independent noise sources, but 
still satisfies Eq. 1261 with the same matrix T. For that 
reason, it can be assumed without loss of generality that 
M < N. The components of the noise in diagonal coor- 
dinates can be written as 



(27) 



with a set of constants computed from the matrix g 
and the coordinate transform between position-velocity 
coordinates and diagonal coordinates. 

With the noise term, Eg. 1271 the time correlation func- 
tion of two stable components l|20(l of the TS trajectory, 
with Rcej < 0, Ree^ < 0, reads 



la trr J-> 



dre £ ' (t-T) 



dae-**" (W am (T)W an (a)) a . (28) 



Using Eq. for t > and with 

^jk — ^ ^jn^kn j 
n 

this can be evaluated to 



(29) 



Similarly, 



and 



if Reej < 0, Ree k < 0. (30a) 



a €j + Cfc 

if Reej > 0, Ree fc > 0, (30b) 



/a e 7 + e/c 



if Reej < 0, Ree fc > 0, (30c) 



4j(*)4fe(°)> = if Reej > 0, Ree fc < 0. (30d) 
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FIG. 3: The instance of the TS trajectory shown in Fig. [5] and a transition path under the influence of the same noise. 
Trajectories are projected onto (a) the reactive degree of freedom, (b) the transverse degree of freedom, and (c) configuration 
space. At the lower end of each column is a projection into the corresponding phase-space plane, disregarding time. The 
instantaneous positions of the moving coordinate axes (dotted) and, in (a), of the invariant manifolds (dashed) are shown at 
t = 0, at t = 22, and at the unique reaction time t rcac t = 8.500 where the TS is crossed. The solid line at i rea ct indicates the 
moving TS The curves in the top faces of each panel illustrates the dynamics of the transition path in relative coordinates Aq, 
Av (not to scale in (a), for graphical reasons). 



All these correlation functions decay exponentially for 
long times. Note that that the correlation in Ea. l30dl of 
an unstable component at a later time t > with a stable 
component at time zero vanishes identically because the 
former depends on the fluctuating force after time t, the 
latter on the force before time zero. For the same reason, 
the correlation in Eq. l3L)cl of a stable component at t > 
with an unstable component at time zero vanishes as 
t -► 0. 

The correlation functions in Eq. I3UI can be evaluated 
explicitly for one-dimensional dynamics. In this case, the 
eigenvalues are given by Eq. 1151 and the matrix A can be 
diagonalized by 



with 

M=( l M and M- = ^-f £ " ~A . 

In one dimension, the fluctuating force in Eq. 1241 can 
have only one independent component with a scalar cou- 
pling strength g, and using Eq. EH the diagonal coupling 
strength reduces to 



Ki . /ly ^Cg — . (33) 



We will in the following only discuss the correlation 
A — M ( €s ^ \ M^ 1 (31) functions in Eq. [20] for vanishing time shift t — 0. The 

I e u J ^ ' covariances of the stable and unstable components of one- 
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1 r- 

0.9 
0.8 - 




1 2 3 4 5 

Y /(2<o b ) 



FIG. 4: Variances il' J > II of stable (solid curve) and unstable 
(dashed curve) components the TS trajectory z\(t) in one 
dimension. 




1 2 3 4 5 

Y /(2o) b ) 



FIG. 5: Variance ffifip of the position coordinate of the TS 
trajectory zj(t) in one dimension. 



dimensional dynamics reduce to 

K 2 

la Of ' 



(44) a 



-2e s 



. 



(34a) 

(34b) 
(34c) 



The variances, Eq.[«2l of the stable and unstable compo- 
nents are shown in Fig. 01 as a function of the damping 
constant 7, which is proportional to the noise strength. 
Not surprisingly, both variances vanish in the absence of 
noise (7 = 0). In the limit of strong noise, the variance 
of the stable component tends to zero as I/7 2 , whereas 
the variance of the unstable component approaches its 
maximum value k^T '/uj^. 

Using the coordinate transformation in Eq. 1321 the di- 
agonal covariances in Eq. 1341 can be converted into the 
covariances of the position q\ and velocity v\ of the TS 
trajectory. They are 



2 

2 v^2 



(4 2 ) a = ^ 



with 



(35a) 
(35b) 
(35c) 



(36) 



Somewhat surprisingly, the position and velocity of the 
TS trajectory turn out to be statistically independent. 
The variance of the position coordinate is shown in Fig. [5] 
as a function of 7. Similar to the behavior of the unstable 
component, it grows with increasing noise strength from 
zero to a maximum value of k^T /uj 2 ^. 

These geometrical observations lead to a novel inter- 
pretation of the friction dependence of reaction rates: In 



a well-defined reactive system there is a potential well 
with a stationary population of rcactants on one side of 
the barrier. Only the fraction of the population that lies 
beyond the stable manifold of the moving reference point, 
which is very small for small friction, will be reactive. As 
the friction increases, the reference point moves further 
down into the reactant well so that a larger portion of 
the population can at times become reactive. That ef- 
fect, however, does not grow indefinitely as the friction 
gets large. Instead, it saturates, so that an opposing 
effect takes over: As was noted in section HV Al the in- 
variant manifolds of the relative dynamics rotate toward 
the coordinate axes for large friction, so that the reactive 
portion of the relative phase space shrinks to zero. This 
effect eventually leads to turnover and to a decreasing 
overall reaction rate at very large friction. 



D. Stochastic separatrices 

The previous construction provides a separation be- 
tween the dynamics of the barrier crossing, as represented 
by the relative dynamics, and its statistical properties, 
which are embodied in the stochastic TS Trajectory. The 
combination of these two aspects allows one to obtain a 
detailed picture of the barrier crossing. As an example, 
we present the calculation of the stochastic separatrix 
near the barrier top, i.e. the collection of all points in 
phase space for which the reaction probability is 50%. 37 

To this end, pick a point z in phase space as the initial 
point, at t = 0, of a stochastic trajectory. This point can 
be characterized equivalently either by its coordinates 
<f(0), v(0) or by the diagonal coordinates z s (0),z u (0) in 
the stable and unstable directions. The latter represen- 
tation together consist of the unstable degree of freedom 
— illustrated in figure ^a), — and an arbitrary number 
of stable transverse coordinates. The fate of a trajec- 
tory in the remote future is determined by the relative 
coordinate Az u (0) = z u (0) — z| u (0): A trajectory will 
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evolve into products as t — > oo if Az u (0) > 0, i.e. if 
z iu(0) < z u (0). Because ^| u (0) is Gaussian distributed 
with zero mean and a variance a\ given by (|30b|l or, in 
one degree of freedom, by (|34a[l . this probability is easily 
computed to be 



P + = I erfc 



(37) 



in terms of the complementary error function 50 
2 f°° 

erfc(x) = —= exp(-t 2 ) dt . (38) 
V n Jx 

Similarly, a trajectory will come from the product side in 
the remote past t — > — oo if z|, s (0) < z s (0), the probabil- 
ity of which is 



P =ierfcf-i^I, 



(39) 



with the variance of given by Ea. I30al or Ea. 134 

Equations 1371 and 1391 provide a means to identify the 
stochastic separatrices, where the probabilities for a tra- 
jectory to evolve into reactants or products are both 
equal to 50%. The stochastic separatrix between trajec- 
tories that evolve into reactants or products in the future 
is the subspace z u (0) = according to Eq. [37| whereas 
the separatrix between reactants and products in the past 
is 2 S (0) = 0. The present considerations thus reproduce, 
in a more quantitative way, the result stated in Ref. [37] 
and at the same time generalize it to systems with ar- 
bitrarily many degrees of freedom. Note, however, that 
the stochastic separatrices separate reactive from non- 
reactive trajectories only in a probabilistic sense, whereas 
the moving separatrices introduces in Sec. II VAl allow one 
to determine the fate of a trajectory with certainty. 



E. Sampling the TS trajectory 

The timc-corrclation functions in Eq. 1301 completely 
specify the probability distribution of the TS trajectory, 
i.e. the joint distribution of all components of the tra- 
jectory at all different times. To use the TS trajectory 
in numerical work of any kind, it is essential to obtain 
random samples z\{t) of the trajectory from this distri- 
bution. The obvious way to do so, namely to generate 
an instance of the fluctuating force £ Q (t) and then com- 
pute the integrals in Eas. 1201 and 1231 would quickly turn 
out to be prohibitively laborious, so that a more efficient 
algorithm must be devised. 

The problem at hand can be states more precisely as 
follows: To sample a representative of the TS trajectory 
z\(t) with L equidistant time steps over a time interval 
T = L At, one has to compute L + 1 random vectors 
4(0), zl{At), . . . , i|(T = L At) from a Gaussian distri- 
bution with zero mean and covariances given by Eq. 1301 
Phrased this way, the problem is to sample a correlated 



Gaussian random variable with 2N(L + 1) components. 
This is easy to do if only a single sample point z|(0) on 
the TS trajectory is required, but if the number of sample 
points is large, it is still a formidable task. 

In the case of white noise, this task is greatly simplified 
by the recurrence relation 



4(t + *) = e e " (4-(i) + A QJ (M)) 



(40) 



where 



A aj (t,s) 



l+s 



^■(r)dr , 



(41) 



that can easily be derived from either Eq. [2U] or Eq. [23 
and holds for both the stable and the unstable compo- 
nents of z\{t). As the fluctuating force £ a j(t), the incre- 
ment functions A a j(t, s) are Gaussian random variables 
with zero mean. Their covariances can be computed from 
Eqs. 03 and 03 to be 



(A aj (t,s) A.* (*,«)>. 



K 



jk 



(42) 



where Kjk is given by Eq. [25] 

Because for white noise the values of the fluctuating 
force at different times are independent, the increment 
functions A a j(ti, si) and A a k{t2, S2) are stochastically 
independent if the time intervals (t%, t\ + si) and (£2, £2 + 
S2) do not overlap, i.e. if t\ + s\ < t% or ti + S2 < 
t\. Further, A a j(t, s) is independent of the component 

z ak( T ) 01 tne TS trajectory whenever either k denotes a 
stable direction with Ree^ < and r < t or k denotes 
an unstable direction with Reefe > and r > t + s. 

These observations lead to the following algorithm for 
sampling the TS trajectory: Generate random values for 
the stable components z a ^ (0) at time zero using the corre- 
lation function, Eq. l30al Generate values for the unstable 
component z^ (T) at time T. They are independent of 
the stable components computed before. Typically, there 
will be only one stable component. Generate increment 
vectors A a (l At, At) for I = 0, 1, . . . , L — 1 using the co- 
variances, Ea. 1421 Increments for different I are indepen- 
dent. Using these increments and the recurrence relation 
in Eq. ^UJ compute the stable components z a j (I At) for 
/ = 1, 2, . . . , L. This recursion is numerically stable be- 
cause the exponential factor in Eq. 0UJ has an absolute 
value less than one. The unstable components z a j(l At) 
are computed using the same recurrence relation, Eq. 1401 
in the opposite order, starting from I = L — 1,...,0. 
Again, the recursion is numerically stable. This algo- 
rithm allows one to sample an arbitrary stretch of the 
TS trajectory without ever having to compute correlated 
random variables of a dimension larger than the phase 
space dimension 2A^. 
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V. THE TRANSITION STATE IN COLORED 
NOISE 

It might seem at first sight that in the derivation of 
the TS trajectory in section llVBI we have never used the 
assumption that the fluctuating force represents 
white noise and that, therefore, the results in Eqs. [201 
and [231 for 5* (t) are equally valid if the noise is col- 
ored. However, if the values of the fluctuating force 
£a(t) a t different times are statistically correlated, the 
fluctuation-dissipation theorem demands the occurrence 
of memory friction effects, and the Langevin equation, 
Eq. 01 must be replaced with the generalized Langevin 
equation 38 i 39 i 41 i 42 i 43 i 44 i 45 i 51 

&(«) = □&(*)- f dTT{t-T%{T) +£,(<) , (43) 

J — oo 

where the symmetric-matrix valued function T(t) is re- 
lated to the fluctuating force by 

(Ut)£(s)) a = k B Tr(\t~s\) . (44) 

Defining 

r(t) = if t < 0, (45) 

the upper limit of integration in Ea . 1431 can be extended 
to infinity. 

Because the friction force in Eq. 031 depends on the 
entire prehistory of the trajectory rather than only the 
present velocity, Eq. 221 cannot be solved as an initial 
value problem. Instead, the entire trajectory up to the 
present time must be specified before its further evolu- 
tion can be determined. Although this freedom seems 
to render the phase space of Ea . 1431 infinite-dimensional, 
we will show below that in many cases the phase space 
dimension is actually finite because only those initial tra- 
jectories can be prescribed that themselves satisfy Ea. 14*51 

To solve the GLE, Eq.021 we use the bilateral Laplace 
transform, 

/oo 
dte-*q(t), (46) 
-OO 

and its inverse^ 

i pioo 

q(t) = 7T- dee^e). (47) 

27 ™ J -too 

It transforms Eq. 021 into 

(e 2 + ef(e)-fi)«f a (e)=|;(e). (48) 

Using Eq. 021 the bilateral Laplace transform T of the 
friction kernel coincides with the customary one-sided 
Laplace transform and is holomorphic on a right half 
plane of the complex e-plane that, for bounded T(t), in- 
cludes the half plane Ree > 0. We further assume r(e) 
to be meromorphic on the entire complex plane. 



A. Relative dynamics 

The relative coordinate Aq(t) , which is again given by 
Eq. |21 satisfies the noiseless version of Eq. 021 vrz - 

(e 2 + ef (e) - o) A<f(e) = . (49) 

The general solution of Eq. 021 reads 

A|te)=5>^5(e- ej ), (50) 

3 

corresponding to 

Aq(i) = ]Tw^, (51) 

3 

where Cj are arbitrary constants and tj , Vj denote a com- 
plete linear independent set of solutions of the nonlinear 
eigenvalue equation 

(e] + e j T(e j )-n)v j =Q, t£ = 1 . (52) 

Because the friction kernel T(t) is real, its Laplace trans- 
form satisfies T(e*) = r(e)*, where the asterisk denotes 
complex conjugation, so that the eigenvalues and eigen- 
vectors of Eq. 1521 are either real or occur in complex con- 
jugate pairs. In the latter case, the constants Cj must 
also be chosen in conjugate pairs to make Aq(t) real. 
The solution, Eq. [31 of the relative dynamics is exactly 
analogous to the solution, Eq. [n] in the case of white 
noise, except that the eigenvalues and eigenvectors are 
computed differently. The construction of invariant man- 
ifolds and a no-recrossing surface that we have given in 
section HV Al therefore carries over unchanged to colored 
noise. A trajectory is specified by prescribing the con- 
stants Cj. In other words, the dimension of the phase 
space of the generalized Langevin equation, Eq.021 is the 
number of solutions of the nonlinear eigenvalue Eq. la"""! 
That number is typically larger than the phase space di- 
mension 2N of the memoryless Langevin equation, Eq. 01 
but still finite if the friction kernel decays exponentially. 
The phase space is described implicitly by Eq. IST1 



B. The Transition State Trajectory 

A particular solution of Eq. 021 in the presence of the 
fluctuating force can be found by solving for q(e) using 
the inverse Laplace transform l|47|) : 

1 r io ° / . 

&(t) = W~ / d€eet I 6 ' + - ") ^( £ ) ' 

(53) 



12 



With the definition of £(e) and after the order of integra- 
tion is interchanged, Eg. 1531 takes the form 



si (+\ - 





\— r 


/ dr 


> — oc 


2m J_ 



e 2 + eT(e) - n 



Ur) ■ (54) 



The integral over e can be evaluated using the calculus 
of residues if the contour of integration is closed through 
the left half plane if r < t and through the right half 
plane if r > i. The poles of the integrand are precisely 
the eigenvalues £j, where the matrix of Eq. 1521 is not 
invertible. The residue of the matrix inverse in (|54H at 
e = €j is denoted by fij, so that 



with 



(55) 



fij / dre-^ T ^ a {t + T) 



(56) 



if Re e, < and 



= -fij dre-^^t + T) (57) 
Jo 

if Reej > 0. These solutions are close analogues of 
Eos.l2Ulandl23lin the case of white noise, with the residue 
factors fij replacing the projection onto the eigenspace j. 
The components ■ can easily be verified to satisfy the 
equation of motion 



(58) 



that is expected for a component in an eigenspace of 
the noiseless dynamics with eigenvalue tj. Again, the 
residues fij play the role of the projection onto the 
eigenspace j. 

A more explicit form of the residues fij can be found 
as follows: First, by observing that 



A(e) = e 2 + eT(e) - fi 



(59) 



is a symmetric matrix, we immediately conclude that fij 



(—fij) is also symmetric. For any fixed e, the matrix 



A(e) possesses a complete set of eigenvalues r]i(e) and 
orthonormal eigenvectors Wi(e): 



A(e)wi(e) = r)i(e)wi(e) , 
Wi{e) ■ Wj(e) = 5 i:j . 



(60) 



In terms of its eigensystem, the matrix A(e) can be writ- 
ten as 



A( £ )=^ry 4 ( e M( e K T (e), 



so that 



, »7i(e) 



(61) 



(62) 



When e is equal to a solution tj of the nonlinear eigen- 
value problem in Eq. 1521 one of the eigenvalues r\i van- 
ishes, say 771 (cj) = 0, and w\(ej) = Vj. If all the eigenval- 
ues are simple, only the term i = 1 in Eq. 1621 has a pole 
at e = e, , and its residue is 



fij = KjVjVj 



where 



(63) 



(64) 



and the prime denotes the derivative with respect to e. 
The residue fij is thus proportional to the geometric pro- 
jection onto the eigenvector Vj, with the factor Kj describ- 
ing the strength of the coupling of the noise to the par- 
ticular eigenmode. Using first-order perturbation theory, 
this factor is found to be 



(65) 



2e i +t£(f(e i )+e i r'(e i ))«/ 



The representation, Eq. ESI of the residues leads to a 
more explicit form of the components q^j (t) : 



<f aj (*) = / dr e-^vj -Ut + t) (66) 

if Ree^ < 0, and a similar expression for Ree^ > 0. Note 
that q^j (t) is a scalar multiple of the eigenvector Vj . 

The eigenvectors of the linear eigenvalue problem in 
Eq. 1201 satisfy the completeness condition 



(67) 



where I is the NxN unit matrix. Because the eigenvalue 
problem in Eq. 1521 is nonlinear, the analogous relation for 
the Vj does not hold. Instead, if T(e) is meromorphic at 
infinity, the residues fij satisfy the completeness condi- 
tions 



/ , t*3 — 

3 3 

X/ X! 'V''/ 1 

3 3 



(68a) 
(68b) 
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To see this, note that 

E^ = ^f/ e ( e2+ef ( £ )- n )" 1 



(69) 



where C is a contour in the complex e-plane that encircles 
all eigenvalues ej in the clockwise direction. With 9 = 
1/e, Ea. 1551 reads 

Ett = ^ f i/c d0(l + 6 f (1/9) - 9 2 n) ~* , (70) 

where the contour l/C is again oriented clockwise. By 
construction, it encircles none of the inverse eigenvalues 
1/ej. As the integrand is furthermore regular at 6 = 
0, the integral is zero, thus proving Eq. I68al To prove 
Eq. I68bl a similar computation leads to 



l = i i deeL 2 + et{e)-Sl) 



1 



d9 

2tt« Ii/c 
= 1, 



I + 9T (l/6») -9 2 n 



r 



(71) 



where the last equality holds because the integrand in 
Eq-E]has a simple pole at 9 = 0. 



C. Statistical properties of the TS trajectory 

As in the case of white noise, the components of the 
TS trajectory, Eos. 1551 and 1571 are Gaussian distributed 
with zero mean if the noise £ a (t) is. To specify their 
distribution completely, it remains only to compute their 
time correlation functions. For t > and in the case of 
two stable components with Reej < 0, Ree^ < 0, they 



(£.(r)6,(*)) a Mfc 
= k B T ii j I jk {t) n k (72) 

with the matrix-valued correlation integral 
I jk (t)= [ dr j dae^-^e-^TQT-al) . (73) 

J — OO J — OO 

Recall that fi k is symmetric, and note that the time cor- 
relation matrix in Eci. l72l is a scalar multiple of the matrix 
Vjv£ if Eq. 1551 holds. 

To evaluate the correlation integral in Eq. [75j we split 
the range of the r integration into the intervals (— oo, 0) 
and (0, t). The latter part can then be expressed through 
the function 

ft />oc 

C(Ci,C2,i)= dr daeM^e-^Tir + a) . (74) 
Jo Jo 



The former part, apart from a factor e £jt , is 

/o j-Q 
dr dae-^ T e~ eka T(\T-a\) , (75) 
-oo J —oo 

which in terms of x — r — a an d tp = t + a, can be 
rewritten as 



,(«fc- e i)x/2 



xl 



l,k(0) = \ J_JxT(\x\)e 



dip e 



+ / d X e-^T( X ) 



+ r(- efc ; 



The correlation integral in Eq. [75J thus finally reads 

f,e } t 



(76) 



-(ej +e k ) 



r(-ei) + r(-e fc ) + C(e„-e fc ,i) 



if Ree, < 0, Ree fc < 



(77a) 



The correlation functions involving unstable components 
take the same general form of Eq. 1721 with 



I jk (t) = (ffe)+f(e fe )) +C(-e k ,e j ,t) 

6j + e k V / 



£j + £fc 



if Ree,- > 0, Ree fc > 



t(e k )+t{-e k ) 



(77b) 



ej + e k 



T(-ej) +T(-e k ) ) + C{e j -e k ,t) 



if Ree,- < 0, Ree fc > 



(77c) 



Ijk(t) = —— r(e,) - r(-e fe ) 1 +C(e h -e k ,t) 
£j + £fe ^ ' 

if Reej > 0, Ree fc < . (77d) 

All correlation integrals consist of exponentially decay- 
ing contributions and a correction term C. The latter can 
be written as 

C(d,C2,*)= [ dre^-^t^) , (78) 
Jo 

where T T denotes the Laplace transform of the time- 
shifted function T(t + r). Unless T(e) is growing at least 
exponentially as |e| — * oo and Ree < 0, an explicit ex- 
pression for the correction term C can be found. To 
this end, use the definition of IV and the inverse Laplace 
transform, Eq. 07| to write, for Ree > 0, 



fOO -i />ZOO 

f T (e)=/ doe-™—\ d V e^+^f( V ) 

Jo 2 ™ J-too 

1 ^drje^M 



2ni 



e-rj 



(79) 
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If it is now permissible to close the integration contour 
through the left half plane, Eg. 1791 yields 



1 „E iT 



e - Ei 



(80) 



where are the poles of f(e) in the left half plane 
(which typically are all poles) and vi are the correspond- 
ing residues. The correction term C thus reads 



C(Ci,C 2 ,*) = E' 



{Ei-SMEi-b) 
and is also exponentially decaying. 

D. Reconstructing the eigenvalues 

The dynamics of the GLE, Eq. 021 is characterized 
by the eigenvalues ej. If the TS trajectory does in- 
deed capture the essential aspects of the dynamics, it 
should enable us to identify the eigenvalues. In particu- 
lar, we should be able to find the Grote-Hynes reaction 



frequenc; 



,34,51 



viz. the unique eigenvalue with positive 



real part, that approximately describes the influence of 
colored noise on the reaction rate. That this information 
can actually be obtained from the TS trajectory can be 
seen most easily by examining the statistical properties 
of its Laplace transform rather than the time-correlation 
functions of the TS trajectory directly. According to 
Eq.EH 



I(e)=(e 2 + ef(e)-n) ^(e) 



(82) 



This Laplace transform has poles at the eigenvalues e,-. 
In addition, there are the random poles of the Laplace 

transform £ a (e) of the noise. Thus, the eigenvalues of 
the noiseless dynamics cannot be determined from the 
Laplace transform of a single instance of the TS trajec- 
tory. If, however, the Laplace transform of several in- 
stances of q| is given, the eigenvalues can be discerned 
with certainty because they are fixed whereas the poles 

of £a(e) are random. 

The randomness of the poles is entirely absent from 
the correlation function of q^, , which is 



x (V 2 + e'r(e') - fi 



(83) 



If the Laplace transform £ a (e) is replaced with its defi- 
nition in Eq. 1461 the remaining average turns out to be 



£*(e)£(e') 



oo 

dr 
Jo 



T(\r-a\). 

(84) 



In the notation of Sec. IV Gl Eq. E] is the correlation 
integral I ee i (0) between two unstable components. It has 
already been evaluated in Ea. l77bl to be 



6(e)£V) 



r(e) + r( e ') 



a e + e 



(85) 



Therefore, 



-i r( e ) + r(e') 



x e' 2 + eT(e') - Q 



(86) 



and in particular, for e = e', 



-1 r(e) 



x (e 2 + et (e) - Q 



(87) 



This correlation function has poles at e — 0, at the poles 
of r(e) and at the eigenvalues of the noiseless dynam- 
ics. Usually, the friction kernel T will be known, so that 
the eigenvalues tj can be determined from Eq. 1871 On 
the other hand, even if the friction kernel is unknown, 
the eigenvalues can be discerned because they are double 
poles of Eq. El whereas the poles of r(e) will typically 
be simple. In particular, the Grote-Hynes frequency^ 4 *^ 
can be identified as the only pole of the correlation func- 
tion in Eq. E| that has a positive real part. 



VI. CONCLUSIONS AND OUTLOOK 

To generalize the formalism of TST to reactive systems 
driven by noise, we have introduced a time-dependent 
dividing surface that is stochastically moving in phase 
space so that it is crossed once and only once by each 
transition path. This construction requires the identifi- 
cation of a privileged stochastic trajectory that remains 
in the vicinity of the barrier for all times. This TS Tra- 
jectory serves as the origin of a stochastically moving 
coordinate system that carries the no-recrossing TS di- 
viding surface for the noiseless dynamics. The moving 
no-recrossing dividing surface thus obtained describes a 
reaction in the same way as a static dividing surface 
in conventional TST and therefore provides a novel ap- 
proach to the theory of reaction dynamics in a noisy en- 
vironment. The challenge lying ahead is to determine a 
rate formula that makes optimal use of the geometric ob- 
jects associated with the TS trajectory and the moving 
no-recrossing TS dividing surface. 

The further generalization of this approach to all-atom 
representations of solvated chemical reactions is nontriv- 
ial. In the Langevin models, the noise sequence can be 
specified independently of the reaction path or the sol- 
vent modes. For an all-atom solvent, by contrast, this 
is impossible because the noise sequence depends on the 
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trajectory of the chosen coordinates identifying the reac- 
tion: A given solvation-shell structure creates an effec- 
tive cavity that accomodates a particular set of reaction 
geometries more easily than others. Thus a critical ques- 
tion to be pursued in future work is how to identify the 
reaction geometry of the moving TS dividing surface in 
all-atom representations of a solvated chemical reaction. 
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